Understanding deterministic diffusion by correlated random walks 
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Low-dimensional periodic arrays of scatterers with a moving point particle are ideal models for 
studying deterministic diffusion. For such systems the diffusion coefficient is typically an irregular 
function under variation of a control parameter. Here we propose a systematic scheme of how to 
approximate deterministic diffusion coefficients of this kind in terms of correlated random walks. 
We apply this approach to two simple examples which are a one-dimensional map on the line and 
the periodic Lorentz gas. Starting from suitable Green-Kubo formulas we evaluate hierarchies of 
approximations for their parameter-dependent diffusion coefficients. These approximations con- 
verge exactly yielding a straightforward interpretation of the structure of these irregular diffusion 
,_0 ' coeficients in terms of dynamical correlations. 
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i r I. INTRODUCTION 

Q: 

Deterministic diffusion is a prominent topic in the theory of chaotic dynamical systems and in nonequilibrium 
statistical mechanics |l],|2| . To achieve a proper understanding of the mechanism of deterministic diffusion with respect 
to microscopic chaos in the equations of motion, a number of simple model systems was proposed and analyzed. An 
interesting finding was that, for a simple one-dimensional chaotic map on the line, the diffusion coefficient is a fractal 
1 — function of a control parameter ^-^]. Diffusion coefficients exhibiting similar irregularities are also known for more 
complicated sytems such as sawtooth maps jj, standard maps and Harper maps the latter two models 
being closely related to physical systems like the kicked rotor, or a particle moving in a periodic potential under the 
influence of electric and magnetic fields. More recently, related irregularities in the diffusion coefficient were discovered 
in Hamiltonian billiards such as the periodic Lorentz gas Q, and for a particle moving on a corrugated floor under 
the influence of an external field |l(J . That this phenomenon is not specific to diffusion coefficients is suggested by 
results on other transport coefficients such as chemical reaction rates in multibaker map s ]TT[ | , electrical conductivities 
in the driven periodic Lorentz gas fl^| , and magnetoresistances in antidot lattices [fl3[ which, again, exhibit strongly 
irregular behavior under parameter variation. 

In this paper we wish to contribute to the further analysis and understanding of irregular transport coefficients 
in simple model systems. We use certain forms of the Green-Kubo formula as starting points for expansions of the 
diffusion coefficient in terms of correlated random walks. In Section 2 we exemplify this approach for the piecewise 
linear map on the line studied in Refs. |3] 0. In Section 3 we apply a suitably adapted version of it to diffusion 
. in the periodic Lorentz gas. In both cases our approximations converge quickly to the precise diffusion coefficient 
as obtained from other numerical methods. We argue that for both models this approximation scheme provides a 
physical explanation for the irregular structure of these diffusion coefficients in terms of dynamical correlations, or 
memory effects. In Section 4 we summarize our results and relate them to previous work on irregular and fractal 
transport coefficients. 
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II. DETERMINISTIC DIFFUSION IN ONE-DIMENSIONAL MAPS ON THE LINE 



Probably the most simple models exhibiting deterministic diffusion are one-dimensional maps defined by the equa- 
tion of motion 

x n+ i = M a (x n ) (I) 

where a £ M is a control parameter and x n is the position of a point particle at discrete time n. M a (x) is continued 
periodically beyond the interval [0, 1) onto the real line by a lift of degree one, M a (x + 1) = M a (x) + 1. We assume 
that M a {x) is anti-symmetric with respect to x = 0, M a (x) — —M a {—x). Let 
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m a (x) := M a (x) mod 1 (2) 

be the reduced map related to M a (x). This map governs the dynamics on the unit interval according to x n = m"(a;). 
Let p n {x) be the probability density on the unit interval of an ensemble of moving particles starting at initial conditions 
x = xq . This density evolves according to the Frobenius-Perron continuity equation 

p n +i(x)= dy pn(y) 5(x - m a (y)) . (3) 
Jo 

Here we are interested in the deterministic diffusion coefficient defined by the Einstein formula 

D(a) = lim < (x n ~ x ) 2 > l{2n) (4) 

n — too 

with x n governed by M a (x), where the brackets denote an average over the initial values Xq with respect to the 
invariant probability density on the unit interval, (...) := f Q dx p*(x).... Eq.(Q) can be transformed onto the 
Green-Kubo formula for maps [E]]lJl4,Eq] , 



D(a) = 5<jV)>+f>X*o)j(a*)) , (5) 
fc=i 

where the jump velocity 

j(x n ) := [x n+ x] - [x n ] , (6) 

[x] being the largest integer less than x, takes only integer values and denotes how many unit intervals a particle has 
traversed after one iteration starting at x n . The map we study as an example is defined by 

M a (x) = ( ^ aX ' \ <X ~M , (7) 
v ; |_ ax + 1 - a , i < x < 1 J w 

where the uniform slope a > 2 serves as a control parameter. The Lyapunov exponent of this map is given by 
A(a) = In a implying that the dynamics is chaotic. A sketch of this map is shown in Fig. ^ (a). Fig. y furthermore 
depicts the parameter-dependent diffusion coefficient of this map in the regime of 2 < a < 4 as calculated in Refs. 
[]3|-|5| by means of a numerical implementation of analytical methods. In these references numerical evidence was 
provided that the irregularities on a fine scale are reminiscent of an underlying fractal structure. 

To analyze this fractal diffusion coefficient, Eq.(^) forms a suitable starting point because it distinguishes between 
two crucial contributions of the dynamics to the diffusion process, which are (1) the motion of the particle on the 
unit interval x n mod 1 generating the invariant density Pa(x), and (2) the integer jumps from one unit interval to 
another one as related to j(x n ). In fact, in Refs. P,|Tl|] it was shown that both parts are independent sources of 
fractality for the diffusion coefficient. However, in these references Eq. (||) was only discussed in the limit of infinite 
time. In this paper we suggest to look at the contributions of the single terms in this expansion, and to analyze 
how the Green-Kubo formula approaches the exact diffusion coefficient step by step. We do this by systematically 
building up hierarchies of approximate diffusion coefficients. These approximations should be defined such that the 
different dynamical contributions to the diffusion process are properly filtered out. Another issue is how to evaluate 
the single terms of the Green-Kubo expansion on an analytical basis. In this section we show how to do this for the 
one-dimensional model introduced above, in the next section we apply the same idea to billiards such as the periodic 
Lorentz gas. 

We start by looking at the first term in Eq.([|). For a < 4 the absolute value of the jump velocity j(x n ) is either 
zero or one. Assuming that p* a {x) ~ I for a — > 2 and cutting off all higher order-terms in Eq. (|^), the first term leads 
to the well-known random walk approximation of the diffusion coefficient P, |ll] , |l6| -|l8| 

D (a) = ±J o dxf a (x) (8) 

which in case of the map Eq. (R) reads 

D (o) = (o-2)/(2o) . (9) 

This solution is asymptotically correct in the limit of a — > 2. More generally speaking, the reduction of the Green-Kubo 
formula to the first term only is an exact solution for arbitrary parameter values only if all higher-order contributions 
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from the velocity autocorrelation function C(n) :=< j(xo)j(x„) > are strictly zero. This is only true for systems of 
Bernoulli type [Q. Conversely, the series expansion in form of Eq. (|^) systematically gives access to higher-order 
corrections by including higher-order correlations, or memory effects. This leads us to the definition of two hierarchies 
of correlated random walk diffusion coefficients: 

(1) Again, we make the approximation that p*(x) ~ 1. We then define 



D J n (a) :=- dx f{x) + V / dx j(x)j(x k ) , n > , (10) 



with -Do (a) given by Eq. (||). Obviously, this series cannot converge to the exact D{a). 
(2) By using the exact invariant density in the averages of Eq.(H) we define 



£><»:=- / dx p* a (x)f a (x) + J2 / 
1 J ° fc=i J o 



dx P* a (x)ja(x)j a (x k ) , n > , (11) 



here with -Dq(o) = | Jq 1 p* a {x)j1(x) 1 which of course must converge exactly. The approximations D 3 n (a) and DP(a) 
may be understood as time-dependent diffusion coefficients according to the Green-Kubo formula Eq. (||) . According 
to their definitions, D J n (a) enables us to look at contributions coming from j(x n ) only, whereas D^(a) assesses the 
importance of contributions resulting from p* a (x). The rates of convergence of both approximations give an estimate 
of how important higher-order correlations are in the different parameter regions of the diffusion coefficient D(a). 
Similarly to Eq. ([)]), the next order -Di(a) can easily be calculated analytically and reads 

!(o-2)/(2a) , 2<a<l + \/3 and 2 + V2<a<4 > | 
3/2-3/a-2/a 2 , 1 + %/3 < a < 3 > . (12) 

-l/2 + 3/a-2/a 2 , 3 < a < 2 + V2 J 

Further corrections up to order n = 3 were obtained from computer simulations, that is, an ensemble of point 
particles was iterated numerically according to Eq. (Q) . All results are contained in Fig. [j] showing that this hierarchy 
of correlated random walks generates a self-affine structure which resembles, to some extent, the one of the well-known 
Koch curve. Fig. [j] (a) illustrates that this structure forms an important ingredience of the exact diffusion coefficient 
D{a) thus explaining basic features of its fractality. Indeed, a suitable generalization of this approach in the limit of 
time to infinity leads to the formulation of D(a) in terms of fractal generalized Takagi functions [p|,pd]. 

Fig. |l| (b) depicts the results for the series of D£(a) up to order n = 2 as obtained from computer simulations. This 
figure illustrates that there exists a second source for an irregular structure related to the integration over the invariant 
density, as was explained above. In the Green-Kubo formula Eq.(||) both contributions are intimately coupled with 
each other via the integration over phase space. 

We now perform a more detailed analysis to reveal the precise origin of the hierarchy of peaks in Fig. |l|. For this 
purpose we redefine Eq. (|l(]) as 

Di(a)= [ dxj{x)J n {x)-\ ( dxf{x),n>0 , (13) 
Jo z Jo 

where we have introduced the jump velocity function 

n 

J n (x) ;=J2j(zk) ■ (14) 

fc=0 

From Eq. it follows J n (x) — [x n +i], that is, this function gives the integer value of the displacement of a particle 
starting at some initial position x. In Fig. || we depict some representative results for J\{x) under variation of the 
control parameter a. Because of the symmetry of the map we restrict our considerations to < x < 0.5. Eq. (|l3|) 
tells us that the product of this function with j(x) determines the diffusion coefficient Dia). The shaded bar in Fig. |2j 
marks the subinterval in which j(x) = 1, whereas j(x) — otherwise, thus an integration over J\(x) on this subinterval 
yields the respective part of the diffusion coefficient. One can now relate the four diagrams (a) to (d) in Fig. |2| to the 
functional form of D\{a) in Fig. [j] (a) thus understanding where the large peak in D{(a) for 2.732 < a < 3.414 comes 
from: For a < 2.732, J\(x) does not change its structure and the interval where particles escape to other unit intervals 
increases monotonously, therefore D(a) increases smoothly. However, starting from a — 2.732 particles can jump for 
the first time to next nearest neighbors within two time steps, as is visible in J\(x) taking values of 2 for x close to 
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0.5. Consequently, the slope of D(a) increases drastically leading to the first large peak around a = 3. Precisely at 
a = 3, backscattering sets in meaning that particles starting around x = 0.5 jump back to the original unit interval 
within two time steps, as is reminiscent in Jx(x) in form of the region J±(x) = for x close to 0.5. This leads to the 
negative slope in D(a) above a — 3. Finally, particles starting around x — 0.5 jump to nearest neighbor unit intervals 
by staying there during the second time step instead of jumping back. This yields again a monotonously increasing 
D(a) for a > 3.414. Any higher-order peak for D J n (a),n > 1 , follows from analogous arguments. Thus, the source 
of this type of fractality in the diffusion coefficient is clearly identified in terms of the topological instability of the 
function J\{x) under variation of the control parameter a. Indeed, this argument not only quantifies two previous 
heuristic interpretations of the structure of D(a) as outlined in Refs. |^-|e|, it also explains why, on a very fine scale, 
there are still deviations between these results and the precise location of the extrema in D(a), cp. to the "overhang" 
at a = 3 as an example. The obvious reason is that contributions from the invariant density slightly modifying this 
structure are not taken into account. 

Looking at the quantities J n (x) furthermore helps us to learn about the rates of convergence of the approximations 
D J n (a) to D{a) at fixed values of a, as is illustrated in Fig. ||| (a) to (d). Here we have numerically calculated J n (x) 
at a = 3.8 for n = 0, 1, 2, 3. Again, the shaded bar indicates the region where j(x) = 1 enabling J„(x) to contribute 
to the value of the diffusion coefficient according to Eq. (|l3|). In fact, J n (x) may also be interpreted as the scattering 
function of an ensemble of particles starting from the unit interval, since it sensitively measures the final position to 
which a particle moves within n time steps under variation of its initial position x. One can clearly see that, with 
larger n, J n {x) develops more and more discontinuities eventually leading to a highly singular and irregular function 
of x. Integration over further and further refinements of J n (x) determines the convergence of the series of D 3 n {a) to 
a fixed value D^a), cp. to Fig. [l] (a). 

We remark that a suitable integration over the functions J n (x) leads to the definition of fractal so-called generalized 
Takagi functions, which can be calculated in terms of de Rham-type functional recursion relations In a way, 

the integration over jump velocity functions such as the one shown in Fig. |3| (d) is similar to the integration over 
Cantor set structures leading to Devils' staircase- type functions |llj . Our results presented so far thus bridge the gap 
between understanding the coarse functional behavior of D(a) on the basis of simple random walk approximations 
only, and analyzing its full fractal structure in terms of Takagi-like fractal forms, in combination with an integration 
over a complicated non-uniform invariant density. We now show that essentially the same line of argument can be 
successfully applied to more physical dynamical systems such as particle billiards. 



III. DETERMINISTIC DIFFUSION IN BILLIARDS 



The class of two-dimensional billiards we want to discuss here is described as follows: A point particle undergoes 
clastic collisions with obstacles of the same size and shape whose centers are fixed on a triangular lattice. There is no 
external field, thus the equations of motion are defined by the Hamiltonian H — rav 2 /2 supplemented by geometric 
boundary conditions as induced by the scatterers. A standard example is the periodic Lorentz gas for which the 
scatterers consist of hard disks of radius R, see Fig. [| In the following we choose m = 1, v = 1, R = 1, and 

as a control parameter we introduce the smallest inter disk distance w such that the lattice spacing of the disks 
is 2 + w. w is related to the number density n of the disks by n(w) = 2/[V3(2 + w) 2 ]. At close packing w = 
the moving particle is trapped in a single triangular region formed between three disks, see Fig. ||, part (1). For 
< w < Woo — 4/v / 3 — 2 = 0.3094, the particle can move across the entire lattice, but it cannot move collision-free 
for an infinite time. For w > Woo the particle can move ballistic-like in form of arbitrarily far jumps between two 
collisions. 

The diffusion coefficient for this particle billiard can be defined by the two-dimensional equivalent of the Einstein 
formula Eq.(|J) reading 

D(w) = lim < (x(t) - x(0)) 2 > /{At) , (15) 

t — >oo 

where, again, the average is taken over the equilibrium distribution of particles with position coordinates x(£). It can 
be proven that in the regime of < w < w M the parameter-dependent diffusion coefficient D(w) exists p0[ . The full 
parameter-dependence of this function was discussed particularly in Ref. || showing that, on a fine scale, D(w) is 
again an irregular function of the parameter w similarly to the diffusion coefficient of the one-dmensional map D(a) 
as discussed above. Whether D{w) is indeed fractal, or maybe C 1 but not C 2 in contrast to the one-dimensional map 
discussed above, is currently an open question. The main issue we want to focus on in this section are quantitative 
approximations for the full parameter dependence of D(w), and to check for the importance of memory effects. A 
first simple analytical approximation for the diffusion coefficient was derived by Machta and Zwanzig in Ref. pT| . 
This solution was based on the assumption that diffusion can be treated as a Markovian hopping process between the 
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triangular trapping regions indicated in Fig. |]. For random walks on two-dimensional isotropic lattices the diffusion 
coefficient then reads D = £ 2 /(4r), where £ = (2 + w)/y/3 is the smallest distance between two centers of the traps, 
and t -1 is the average rate at which a particle leaves a trap. This rate can be calculated by the fraction of phase 
space volume being available for leaving the trap divided by the total phase space volume of the trap thus leading to 
the Machta-Zwanzig random walk approximation of the diffusion coefficient 

n / ^ w(2 + w) 2 

£>mzH = r k 10 , , 2 - 1 ■ (16) 

Indeed, this approximation is precisely the billiard analogue to the one-dimensional random walk diffusion coefficient 
for maps Eqs.(^),(^|). Similarly, this approximation is asymptotically exact only for w — > 0, as is shown in comparison 
to computer simulation results in Fig. || [ppl). 

In Ref. ||, Eq. ( |l6|) was systematically improved by including higher-order correlations. Two basic approaches were 
presented both starting from the idea of Machta and Zwanzig of looking at diffusion in the Lorentz gas as a hopping 
process on a hexagonal lattice of "traps" with frequency t _1 . This picture was quantified by introducing a simple 
symbolic dynamics for a particle moving from trap to trap as indicated in Fig. ^. Let us follow a long trajectory of 
a particle starting with velocity v parallel to the y-axis, cp. part (1) in Fig. ||. For each visited trap we label the 
entrance through which the particle entered with z, the exit to the left of this entrance with /, and the one to the right 
with r, see part (2) in Fig. 0. Thus, a trajectory in the Lorentz gas is mapped onto words composed of the alphabet 
{z,l,r}. One can now associate transition probabilities to these symbol sequences reading, for a time interval of 2r, 
p(z),p(l),p(r). p{z) corresponds to the probability of backscattering, whereas p(l) = p(r) = (1 — p(z))/2 indicates 
forward scattering. For a time interval of 3r, we have nine symbol sequences each consisting of two symbols leading 
to the probabilities p(zz), p(zl), p{zr) 1 p(lz), p(ll), p(lr), p{rz) 1 p(rl), p{rr), and so on. Beside this hierarchy of 
conditional probabilities defined on a symbolic dynamics there is a different type of probability that a particle leaving 
a trap jumps without any collision directly to the next nearest neighbor trap. As was shown in Ref. @, Eq. @ 
can be corrected by analytically including all these probabilities. Alternatively, lattice gas computer simulations were 
performed by using the probabilities as associated to the symbol sequences. The heuristic corrections of Eq. (|l^) led 
to a satisfactorily explanation of the overall behavior of D(w) on a coarse scale, however, the convergence was not 
exact. The lattice gas simulations, on the other hand, were converging exactly, however, here a proper analytical 
expression for the diffusion coefficient approximations in terms of the associated probabilities was not available. 

In analogy to the procedure as outlined for the one-dimensional map, that is, starting from a suitable Green-Kubo 
formula, we will now define a third approximation scheme which we expect to be generally applicable to diffusion 
in particle billiards. Compared to the two existing approaches mentioned above the advantage of the new method 
is two- fold, namely (1) that by using the set of symbolic probabilities the respective Green-Kubo formula can be 
evaluated according to an analytical scheme, and (2) that the resulting approximations converge exactly to the 
computer simulation results. 

In the Appendix we prove that, starting from the Einstein formula Eq. (|l5|), quite in analogy to the one-dimensional 
case a Green-Kubo formula can be derived which is defined for an ensemble of particles moving on the hexagonal 
lattice of traps depicted in Fig. ||. The result reads 

^H = ^(j 2 ( x o)) + ^E<j(xo)-j(x„)) . (17) 

n=l 

Here j(x n ) defines jumps at the nth time step in terms of the lattice vectors £ Q| a 7 ... associated to the respective symbol 
sequence of the full trajectory on the hexagonal lattice, cp. to Fig. || Let us start with j(x ) = 1/t , I := (0,£)*. 
The next jumps are then defined by j(xi) = £ a /T , a £ {I, r, z}, and so on. The averages indicated in Eq. ( p"7| ) by the 
brackets are calculated by weighting the respective scalar products of lattice vectors with the corresponding conditional 



probablities p(a(3"f . . .). Note that Eq. (17) is the honeycomb lattice an alog ue to the Green-Kubo formula derived 
by Gaspard for the Poincare-Birkhoff map of the periodic Lorentz gas P,p2|. The Poincare-Birkhoff version is very 
efficient for numerical computations, however, according to its construction it fails to reproduce the Machta-Zwanzig 
approximation Eq. (|l6|). Consequently, it does not appear to be very suitable for diffusion coefficient approximations 
of low order. More details will be discussed elsewhere |p3| . We remark that, in terms of using a symbolic dynamics, 
there is also some link between Eq. (|l^) and respective diffusion coefficient formulas obtained from periodic orbit 
theory @,||. 

We now demonstrate how Eq. ( |l7| ) can be used for systematic improvements of the diffusion coefficient on the 
lattice of traps by including dynamical correlations: As in case of one-dimensional maps we start by looking at the 
first term in Eq. ( |l7j ) and cut off all higher order contributions. Obviously, the first term is again the random walk 
expression for the diffusion coefficient on the hexagonal lattice of traps which, by including the respective solution for 
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the jump frequency r 1 , boils down to Eq. (|T^). For calculating higher-order corrections we now define the hierarchy 
of approximations 

D n (w) = L + ±Y,p( a fa...)l.l( a p>y...),n>0 , (18) 



4r 2r 

ap-y... 

with Do(w) given by Eq. ([!(]). To our knowledge yet there is no method available to analytically calculate the 
conditional probabilities p(a(3j . . .). Our following evaluations are therefore based on the data presented in Ref. || 
as obtained from computer simulations. In terms of the formal probabilities it is now easy to calculate the solution 
for the first order approximation at time step 2r to 

D^w) = D (w) + D {w)(l - 3p(z)) . (19) 

For a comparison of this formula with the simulation data D(w) see Fig. [|. We remark that the corresponding solution 
in Ref. || as obtained from a heuristical correction of Eq. ( |l6| ) reads £>i,mz(w) = D (w)3(l — p(z))/2. Indeed, one 
can show that the new formula Eq.([l9]) is closer to D(w) for large w, whereas for small w the previous approximation 
is somewhat better. It is straightforward to calculate the two approximations of next highest order to 



D 2 (w) = Di(w) + D Q {w) {2p(zz) + 4p(lr) - 2p(K) - 4p(lz)) (20) 



and to 



D 3 (w) = D 2 (w) + D (w) \p(llr) + plllz) + p{lrl) + p(lrr) + p(lzl) + p(lzz) + p(rll) + p(rlr) + p{rrl) (21) 
+p(rrz) +p(rzr) +p(rzz) + p(zll) +p(zlz) + p(zrr) +p(zrz) +p(zzl) + p(zzr) 
—2{p{lll) +p(lrz) +p(lzr) +p(rZz) + p(rrr) +p(rzl) + p(zlr) + p(zrl) + p(zzz))] 

All results are shown in Fig. [| demonstrating that the series of approximations defined by Eq. (ffsh converges quickly 
and everywhere to the simulation results. Our new scheme thus eliminates the deficiency of the semi- analytical 
approximation proposed in Ref. |1 that was based on heuristically correcting the Machta-Zwanzig approximation Eq. 
(|l6|). By comparing this new scheme with the lattice gas simulations of the same reference, on the other hand, it 
turns out that the rate of convergence of the lattice gas approach is still a bit better. In any case, all these three 
methods unambiguously demonstrate that for achieving a complete understanding of the density-dependent diffusion 
coefficient in the periodic Lorentz gas it is unavoidable to take higher-order correlations, or the impact of memory 
effects, properly into account. 

We finally remark that a very good low-order approximation for the diffusion coefficient can already be obtained by 
combining Eq. ( |l9| ) with the probability of coUisionless flights p c { (w) mentioned above, that is, by taking into account 
the possibility of next nearest neighbor jumps. The correction of Do(w) as given by Eq. (|l6| ) according to coUisionless 
flights only was already calculated in Ref. |9| and read -D 0jC f(u>) = Dq(w)(1 + 2p c {(w)). Adding now the second term 
of Eq. (|1S|) to this expression by just following the Green-Kubo scheme yields 

£> 1 ,cfM=A>(uO(2 + 2p cf (uO-3p(2)) ■ (22) 

This solution is also depicted in Fig. || and shows that this approximation indeed significantly improves Eq. ( |l9| ) 
for large w yielding a function that is qualitatively and quantitatively very close to D(w). We know of no better 
approximation for D(w) based on information such as p(z) and p c f(w) only. The succesful application of Eq. ( p^ ) 
suggests that coUisionless flights form an important mechanism to understand the full diffusive dynamics of this 
billiard. However, somewhat surprisingly they are not explicitly contained neither in the Green-Kubo expansion Eq. 
( |l8| ) nor in the lattice gas simulations of Ref. || . In both cases exact convergence is achieved by following the hierarchy 
of symbol sequence probabilities only, where coUisionless flights are not apparent. 



IV. SUMMARY AND CONCLUSIONS 



In this paper we suggested a general scheme of how to understand the structure of parameter-dependent determin- 
istic diffusion coefficients in simple model systems. The important point was to find suitable Green-Kubo formulas 
to start with. We used the fact that the class of models we studied here was defined on periodic lattices, and re- 
spectively we discretized the dynamics of the moving particles according to these lattices. As two different examples 
we analyzed a simple one-dimensional map on the line as well as the periodic Lorentz gas. In both cases we recov- 
ered the respective well-known random walk formulas for the diffusion coefficient as zero-order approximations in our 
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Green-Kubo approach. Wc then calculated higher-order terms according to Green-Kubo thus systematically includ- 
ing higher-order correlations. As much as possible this was performed analytically, alternatively in combination with 
data obtained from computer simulations. Our results provided clear evidence that a proper understanding of the 
parameter-dependent diffusion coefficients in both models can only be achieved by taking strong memory effects into 
account. In this respect our research appears to be somewhat related to the findings of long-time tails in the velocity 
autocorrelation function of simple model systems such as Lorentz gases, and to the existence of non-analyticities in 
the density expansion of transport coefficients, which are both consequences of strong correlations in the dynamics 
of the moving particles; see Ref. pq| for a nice review. In case of the one-dimensional map our approach enables 
a detailed understanding of the dynamical mechanism generating the most pronounced irregularities in the fractal 
diffusion coefficient. In case of the periodic Lorentz gas this scheme straightforwardly generalizes the Machta-Zwanzig 
random walk formula in terms of systematic higher-order approximations that converge exactly to the simulation 
results. 

For one-dimensional maps the infinite time limit of this approach, though not the intermediate level as quantita- 
tively discussed here, was already worked out in Refs. j|,[ll| leading to the definition of the diffusion coefficient in 
terms of fractal so-called generalized Takagi functions. For the periodic Lorentz gas an analogous generalization would 
be desirable as well. We furthermore remark that the approximation scheme as presented in this paper was already 
successfully applied (1) to the one-dimensional climbing sine map, where there is a complicated transition scenario 
between normal and anomalous diffusion (2^], and (2) to the so-called flower shape billiard, where the hard disks of the 
Lorentz gas are replaced by obstacles of flower shape |2§| l . In both cases the resulting parameter-dependent diffusion 
coefficients, as far as existent, are much more complicated functions than in the corresponding models discussed above, 
yet our scheme yields systematic explanations of the structure of these functions in terms of strong dynamical corre- 
lations. We thus expect that using Green-Kubo formulas this way provides a general access road to understanding 
deterministic diffusion in low-dimensional periodic arrays of scatterers, possibly also in view of experimental results 
on systems such as antidot lattices, ratchets, and Josephson junctions. 
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APPENDIX A: GREEN-KUBO FORMULA ON THE HEXAGONAL LATTICE 



In this Section we derive Eq. (17) which yields the diffusion coefficient D(w) for the periodic Lorentz gas via a 
generalization of the Machta-Zwanzig picture. That is, we look at diffusion as a higher-order Markov process on a 
hexagonal lattice of traps with lattice spacing & where particles hop with frequency r _1 from trap to trap, cp. to Fig. 
||. The time t is suitably rewritten in terms of the escape time r = (ir/6w)(V3/2(2 + w) 2 — ir) as t = nr , n 6 N. 
Let x ra be the position of the moving particle at time step n. We then write x„ = X„ + x„, where X„ denotes the 
position of the trap on the hexagonal lattice in which the particle is situated at time step n. This vector can be 
expressed by a suitable combination of lattice vectors. For example, one may choose a sum of £ a 3-y... as introduced in 
Section III, see Fig. ^, a more precise definition is not important here. Correspondingly, x n is the distance from the 
nearest trap center to the actual position of the particle in the Wigner-Seitz cell. 

The Einstein formula Eq. (|l5|) then reads 



D(w)= lim <(X„ + Ax„) 2 )/(4nr) 



(Al) 



where Ax n := x n — xo- Multiplying out the nominator we get 



< X 2 + 2X„ • A5 



Ax! > 



(A2) 



According to its definition the last term is bounded, Ax n < const. To the second term we apply the Holder inequality 
M yielding 



< X„ • Ax„ > | < y/< |X„| 2 >V< |Ax„| 2 > < const y/< |X„| 2 > 



(A3) 



Consequently , in the limit of infinite time only the first term in Eq. ( A2) contributes to the positive diffusion coefficient 
D(w) of Eq. (|Al|) leading to 



D(w) = lim < X 2 > /(4nr) 



(A4) 
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Starting from this Einstein formula on the hexagonal lattice it is now straightforward to derive Eq. ([17]). = X^+Y^ 
tells us that essentially there are two one-dimensional parts, thus we are back to the respective derivation of Eq 
in Ref. ffl. To make this precise, let us write 

n-l 

X n = Y,tt*k) . (A5) 

k=0 

which is in analogy to Eq. (|l4|). However, here j(xk) '■= Xk+i — X% is strictly speaking no jump velocity, since it 
only denotes the distance a particle jumps within a time interval r, whereas r = 1 in Eq. (o). Multiplying out the 



nominator of Eq. (A4) in terms of Eq. (A5) leads to an equation for D(w) in form of velocity autocorrelation functions 
C(k,l) :=< j(xk) ■ j(xi) >= J dx dy p*(x,y) j(xk) • j(xi), where p*(x,y) is the equilibrium density of the periodic 
Lorentz gas. Translational invariance implies C(k, I) = C(0,l — k) which is easily shown by substitution combined with 
conservation of probability according to the Frobenius-Perron equation of the billiard. Summing up all contributions 
C(0, k) obtained from the multiplication, doing the same for the component Y n , and putting both results together 
yields Eq. (ET 
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FIG. 1. Diffusion coefficient D(a) for the one- dimensional map shown in the upper left part of (a), where a is the slope of the 
map. The dots are obtained from the method of Refs. [2-4] , the different lines correspond to different levels of approximations 
based on the Green-Kubo formula Eq. (JE]). In (a) the approximations were computed from Eq. (jl(]) assuming a constant 
invariant density, in (b) they are from Eq. ( pj| ) which includes the exact invariant density. Any error bars are smaller than 
visible. All quantities here and in the following figures are without units. 
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X 



X 




FIG. 2. Jump velocity function Ji(x) as denned by Eq. (J1J) which gives the integer value of the displacement of a particle 
starting at some initial position x. Shown are results for different values of the slope a. It is (a) a = 2.7, (b) a = 3.0, (c) 
a = 3.3, (d) a = 3.8. The shaded bar marks the subinterval where the jump velocity j(x n ) defined in Eq. (^) is equal to one. 
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FIG. 3. Jump velocity function J n (x) as denned by Eq. (gj) at the fixed parameter value a = 3.8 for the different number 
of time steps (a) n — 0, (b) n = 1, (c) n = 2, (d) n = 3. Again, the shaded bar marks the subinterval where the jump velocity 
j(x n ) defined in Eq. (^) is equal to one. 
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FIG. 4. Geometry of the periodic Lorentz gas with the gap size w as control parameter. The hatched area related to (1) 
marks a so-called trapping region, the arrow gives the lattice vector connecting the center of this trap to the next one above. In 
(2) three lattice vectors are introduced and labeled with the symbols I, r, z. They indicate the positions where particles move 
along the hexagonal lattice of Wigner-Seitz cells, starting from the trap z, within two time steps of length r. (3) depicts the 
situation after three time steps r with the different lattice vectors associated to symbol sequences of length two. 
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FIG. 5. Diffusion coefficient D(w) for the periodic Lorentz gas as a function of the gap size w as a control parameter. The 
computer simulation results for D(w) axe from Ref. [8], error bars are much smaller than the size of the symbols. The other 
lines correspond to different levels of the approximation Eq. the last approximation (again with symbols) is from Eq. (p2[). 
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